Quantum phase transition using quantum walks in an optical lattice 
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We present an approach using quantum walks (QWs) to redistribute ultracold atoms in an optical 
lattice. Different density profiles of atoms can be obtained by exploiting the controllable properties 
of QWs, such as the variance and the probability distribution in position space using quantum coin 
parameters and engineered noise. The QW evolves the density profile of atoms in a superposition 
of position space resulting in a quadratic speedup of the process of quantum phase transition. We 
also discuss implementation in presently available setups of ultracold atoms in optical lattices. 



I. INTRODUCTION 

Exploiting the aspects of quantum mechanics, such as 
superposition and interference, has lead to the idea of 
quantum walks (QWs), a generalization of classical ran- 
dom walks (CRWs) P, 0, H i, [1] . In the CRW the parti- 
cle moves in the configuration space with a certain prob- 
ability, whereas in the QW the particle moves in a su- 
perposition of the configuration space with a probability 
amplitude. Because of the quantum interference effect, 
the variance in position cr^ of the QW is known to grow 
quadratically with the number of steps N , oc N"^ , 
compared to the linear growth, a'^ (x N, for the CRW. 
This motivated the research in the direction of finding 
quantum algorithms with optimal efficiency using the 
QW [1, 0, S, HI- Experimental implementation of the 
QW has also been reported 0, EH, 111- 

The continuous time quantum walk (CTQW) and the 
discrete time quantum walk (DTQW) are the two wi dely 
studied versions of the QW [H. In the CTQW [11, 
one can directly define the walk on the position space, 
whereas, in the DTQW [15|, it is necessary to introduce 
a quantum coin operation to define the direction in which 
the particle has to move. Due to the coin degree of free- 
dom, the discrete time variant is shown to be more pow- 
erful than the other in some contexts [9)|, and the coin 
parameters (scattering operator) can be varied to con- 
trol the dynamics of the evolution [1, U^. Therefore, in 
this paper we will consider only the discrete time vari- 
ant of the QW. Beyond quantum computation, the QW 
can be used to demonstrate coherent quantum control 
over, for example, atoms, photons or spin chain systems. 
Along with purely quantum dynamics, a small amount 
of engineered noise or environmental effect can enhance 
the properties of the QW [13, [H, [111 ■ As we show later 
this can enlarge the toolbox for controlling the dynamics 
in a physical system. 

This paper describes the use of a QW which consists 
of an iteration of the quantum coin operation and a sub- 
sequent conditional shift operation to control and study 
the dynamics of ultracold bosonic atoms in an optical lat- 
tice. The quantum coin operation evolves the atom wave 
function into the superposition of the internal state of the 
particle and the shift operation spreads the atom wave 
function in a superposition of the position space. In par- 



ticular, we show that the quantum phase transition from 
the Mott insulator (MI) - a regime where no phase coher- 
ence is prevalent - to the superfluid (SF) - a regime with 
long-range phase coherence (20l . [2l| and vice versa can 
be obtained by redistributing the density profile of atoms 
using a QW. The simulation of the quantum phase tran- 
sition using the QW occurs quadratically faster in one 
dimension (ID) compared to, varying the optical lattice 
depth and letting the atom-atom interaction follow the 
CRW behaviour. Enhancing the properties of the QW 
by varying the quantum coin parameters or by addition 
of experimentally engineered noise in small amounts [11 
can be used as an additional tool to evolve the atoms into 
different number density distributions. 

Theoretical studies of the dynamics of atoms in an op- 
tical lattice are done using mean field approaches [Z^ 
and quantum Monte Carlo methods ^23.] . The quantum 
correlation induced between atoms and position space by 
the QW can serve as an alternate method for theoretical 
studies. The use of the QW to study the dynamics of par- 
ticles in ID magnetic systems [24i| , and photonic Mott in- 
sulators |25l| or to observe complex quantum phase tran- 
sitions [26| and quantum annealing [27| would be of both 
theoretical and experimental relevance. 

This paper is organized as follows. In Sec. |TT]we de- 
fine the DTQW and its properties, and in Sec. IIIII the 
phase transition in an optical lattice is discussed. The 
implementation of the QW on atoms in ID MI and SF 
regimes is discussed in Sec. IIVI where we analyze the 
dynamics and present the density profile obtained by ex- 
ploiting some of the properties of the QW. In Sec. fVl we 
discuss the QW with a noisy channel as a tool to control 
the redistribution of atoms. In Sec. I VII we briefly discuss 
a scheme for experimental realization before concluding 
in Sec. Eni 



II. DISCRETE TIME QUANTUM WALK 



To define a ID DTQW we require the 



Hilbert 



space Tic and the position Hilbert space Tip. Tic is 
spanned by the internal state of the particle, |0) and |1), 
and Tip is spanned by the basis states j G Z. The 
total system is then in the space Ti = Tic^Tip. To imple- 
ment the QW, the particle at the origin (j = 0) in state 
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I Jo) is evolved into a superposition of TLc, 
l*in) = ^[|0>+*|l>]»bo) 



(1) 



and then subjected to the conditional shift operation S 
to evolve into a superposition in the position space: 



5 = |o)(o|®Vb--i)(j| + |i)(i| 



|j + i>01 



= |0)(0| ® a- 



Il)(l| 



(2) 



Here a and are annihilation and creation operations 
acting on TLp. The S is followed by the quantum coin 
operation, which in general can be written as an SU (2) 
operator of the form 



cos(6l) e^'^ sin(6l) 



sin(6l) 



cos{0) 



(3) 



to evolve the particle into a superposition in Tic- There- 
fore, each step of the QW is composed of an apphcation of 
operation B and a subsequent operation S to entangle Tic 
and Hp. The process of application ofW = 5(5^, e.c ^ 1) 
is iterated without resorting to intermediate measure- 
ments to realize large number of steps. 
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FIG. 1: Spread of probability distribution for different values 
of 6 using the quantum coin operator Bo.e.o- The distribution 
is wider for (a) = (0, ^,0) than for (b) (0, f , 0) and (b) = 
(0, f§,0), showing the decrease in spread with increase in 6. 
The initial state of the particle \'i'i„) — (|0) (8> |io) and 

the distribution is for 100 steps . 



The variance can be varied by changing the parameter 

K [1 - sm{e)]N^ . (4) 

The effect of the parameter 9 on the distribution is shown 
in Fig. [TJ the variance decreases with the increase in the 
value of 9. The parameters ^ and C introduce asymmetry 
in the position space probability distribution as shown in 
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FIG. 2: Spread of probability distribution for different values 
of 5, 9, <^ using quantum coin operator B^^b^q. Biasing the 
walk using ^ shifts the distribution to the left: (a) = (-1, -I, 0) 
and (c) — -I, 0). Biasing the walk using shifts it to the 

The initial state 



right: (b) = (0,f,f: 
of the particle \^in) - 
is for 100 steps . 



and (d) = (0, 

-- (|0)+i|l))® 
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Fig. [2] and their effect on the variance is negligible [la] ■ 
The simplest coin operation is the Hadamard operator 
H = i3o,45°,o- In this paper. Sec. (|IV|) . we consider the 
QW with i?o,e°,o to redistribute atoms in an optical lat- 
tice. The QW is sensitive to noise [13, El, [3 but in Sec. 
|V]we consider small amounts of three physically relevant 
models of noise that can be experimentally induced: a 
bit-flip channel, a phase- flip channel and an amplitude- 
damping channel [l8|, to act as an enhanced toolbox to 
control the dynamics of the atoms. 



III. PHASE TRANSITION IN OPTICAL 
LATTICE 

A Bose-Einstein condensate (BEC) at low enough tem- 
perature is a SF described b^ a wave function that ex- 
hibits long-range coherence 28] . When the BEC is trans- 
ferred to the lattice potential, the atoms move from one 
lattice site to the next by tunnel coupling. Using a 
second-quantized form, the system can be described by 
the Bose-Hubbard model: 

h = -jY. + E ^^"j- + \uY.^Mo - 1): (5) 

J is the tunneling term characterized by the hopping am- 
plitude, bj and P- are the annihilation and creation oper- 
ators of atoms, Cj denotes the energy offset due to exter- 
nal harmonic confinement of the atoms in the jth lattice 



site, fiA 



I, J — bjbj is the atomic number operator counting 
the number of atoms, and U is the repulsive interaction 
between two atoms in a single lattice site. 
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When J dominates the Hamiltonian, the ground-state 
energy is minimized if the single-particle wave functions 
of all N atoms are spread out over the entire M-lattice 
site. If ej — const, (homogeneous system) then the 
many-body ground state is called the SF state and is 
given by 



where \ V) is a vacuum state. In this state the probability 
distribution for the local occupation rij of atoms on a sin- 
gle lattice site is Poissonian. The state is well described 
by a macroscopic wave function with long-range phase 
coherence throughout the lattice. With increase in the 
ratio U/J, the system reaches a quantum critical point, 
the fluctuations in atom number of a Poisson distribution 
become energetically very costly, and the ground state of 
the system will instead undergo a quantum phase tran- 
sition from the SF state to the MI state, a product of 
local Fock states of n atoms in each lattice site is given 
by MM, 

M 

|^M7)j=ocxn(^]ri^)- (7) 



the atom and position correlations overlap resulting in, 

M 
2 

n if3j-N\0)®\j-N) + 

3 = -^ 

P,-(N+i) |0) ® |j -{N + \)) + + /3,+w|0) ® \j + N) 

+7.-Jv|l) ®\]-N) + +7,+jv)|l> ® I J + A^>IIO) 

This can be written as 

f / 3+N \ 

{wr\^Mi)j=o^ n E [/?.io)+7.ii)]®k) ai) 

j=-M- \x=]-N J 

(3x and 'jx are the probability amplitudes of state |0) 
and |1) at lattice site x, which range from (j — N) 
to {j + N). For a QW of N steps on a particle ini- 
tially at position j = using i?o,e,o as the quantum 
coin, the probability distribution spread over the inter- 
val (— cos(0), iV cos(6')) in position space and shrink 
quickly outside this region [3, HOj. Therefore, after N 
steps the density profile of M atoms initially between 
zb-y will correlate with the position space and spread 
over the lattice site + Ncos{9)). Figure [3] is the 

redistribution of atomic density of 40 atoms initially in 
the MI state when subjected to a QW of different num- 
bers of steps with Hadamard operator H = -Bo,45°.o as 
the quantum coin. 



IV. QUANTUM WALK ON ATOMS IN ID 
MOTT INSULATOR REGIME 



Localized atomic wave functions in the MI regime with 
one atom in each of the M lattice sites are initialized 
into a symmetric superposition of any of the two internal 
trappable state, hyperfine levels, |0) and |1) 29], 



I* 



M/)j=0 « Yl 



|0)+»|1) 
V2 



(8) 



The Tic of each atom is spanned by the two hyperfine 
levels and Hp is spanned by the lattice site. The total 
system is then in the Hilbert space Tim = {HjTicj )®'Hp. 
The unitary shift operation S, Eq. ([2]), on the above 
system will evolve each atom into superposition of the 
neighbor lattice site, establishing the quantum correla- 
tion between the states of the atom and the neighboring 
lattice site. 



S'|^'M7)j=0 OC \\ 



\Q) ® a\j) + ® a^\]) 
V2 



(9) 



To implement the TV number of steps, the process oiW = 
S{B^£^i^ (g) 1) is iterated N times. During the iteration. 
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FIG. 3: Density profile of the evolution of 40 atoms starting 
from MI state in correlation with the position space when 
subjected to a QW of different numbers of steps with the 
Hadamard operator Bo,45°,o as the quantum coin. The dis- 
tribution spreads with increase in the number of steps. 



The redistributed atoms would be in either of their in- 
ternal states and this can be retained to study the phase 
transition of the two-state bosonic atoms in an optical 
lattice 3l| or all atoms can be transferred to one of the 
internal states. A technique based on adiabatic passage 
using crafted laser pulses can be used for a nearly com- 
plete transfer of population between two states [32l | , and 
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a popular example of one such technique is stimulated 
Raman adiabatic passage (STIRAP) [SSj . 

Along with the correlation between the atoms and the 
position space, Eq. PH)) also reveals the overlap of the 
probability amplitude of different atoms in the position 
space. That is, during each step of the QW, the ampli- 
tude of the states of each atom overlaps with the ampli- 
tude of the states of the atoms in the neighboring lattice 



site. After the iteration of the N steps > 



M/2 
cos{9] 



, the over- 



lap of all atoms could be seen at the central region of the 
lattice. 

When the number of steps of the QW is equal to the 
number of lattice site {N = M), the overlap of the frac- 
tion of amplitude of all M atoms exists within the lat- 
tice site ± 



Mc 



s(fl) 



state. Beyond ± 



making the region identical to the SF 
^'^°'^^^'> the number of atoms in corre- 
lation with that position space decreases and hence the 
number of overlapping atoms also decrease. To make all 
M atoms to spread between lattice site using the 
usual method of lowering the optical potential depth fol- 
lowing the CRW protocol takes steps. Therefore, we 
can conclude that using QW a long range correlation can 
be induced quadratically faster than by the usual tech- 
nique. 
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FIG. 4: Density profile of the evolution of 40 atoms start- 
ing from SF state in correlation with the position space when 
subjected to a QW of different numbers of steps with the 
Hadamard operator -Bo,45°,o as the quantum coin. The dis- 
tribution spread with increase in the number of steps and is 
almost uniform between the lattice site ±20 when A*' = 25. 
At this stage the optical potential depth can be increased to 
cancel the correlation and obtain the MI state. 



Similarly, Fig. |4]is the redistribution of atomic density 
of 40 atoms initially in the SF states when subjected to 
the QWs of different number of steps using Hadamard op- 
erator i?o,45°,o as the quantum coin. To reach MI state 
the distribution should be uniform over the lattice site 
without correlation. A distribution that is approximately 
uniform within the region ± ^ can be obtained using the 
QW. In Fig. (g]) when N = 25 the distribution is al- 



most uniform between the lattice site ±20 retaining the 
correlation with the position space. The uniformity of 
the distribution can be improved by using the tunnelling 
of atoms between the lattice. Once the distribution is 
uniform, the optical potential depth can be increased to 
cancel the correlation and obtain the MI state. The uni- 
formity of the distribution can also be improved by in- 
troducing a noise channel, as we show later. 

The variance and the probability distribution can be 
controlled using the parameters 9, ^ and C in B^^g^,^. 
Fig. ([51 shows the density distribution obtained by 
implementing the QW on atoms in MI and SF state 
{N — M — 40) with different values of 6 in the coin 
operator -Bo,0,o- The spread is wider for 9 — 30° and de- 
creases with increase in 9. In Fig. ^ the distribution is 
almost uniform for 9 — 60° between the lattice site ±20. 
At this value the optical potential depth can be increased 
to cancel the correlation and obtain the MI state. 
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FIG. 5: Distribution of atoms initially in MI state when sub- 
jected to a QW (A'^ = AI — 40) using different values for 9 in 
the coin operator -Bo.s.o- The spread is wider for 9 = 30° and 
decreases with increasing 9. 



The Hamiltonian of the system in general can be de- 
scribed by the ID Bose-Hubbard model for two-state 
atoms. 



here "f and J, represent the terms for atoms in state |G) and 
state |1) respectively. U is the interaction between atoms 
in state |0) and |1); V|(|) is the interaction between atoms 
in same state, and dn are the left and right displace- 
ment terms induced during each step of the QW. The 
Hamiltonian is evolved with coin toss operation in a reg- 
ular interval of time t, the time required to move the atom 
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FIG. 6: Distribution of atoms initially in SF state when sub- 
jected to a QW (A'^ = M — 40) using different values for 6 in 
the coin operator Bo.e.o- The distribution is almost uniform 
for 9 — 60° between the lattice site ±20. At this value the op- 
tical potential depth can be increased to cancel the correlation 
and obtain the MI state. 



FIG. 7: Atoms in MI state after implementing the QW 
(A'^ = M = 40,9 = 45°) with noise channel. With increased 
phase damping from pp — 0.02 to pp =0.1 the distribution 
(nj) at the central region gets closure to Gaussian. Ampli- 
tude damping of state |1) introduces asymmetry to the dis- 
tribution, Pa ~ 0.2. 



to the neighboring site. The optical lattice can be dynam- 
ically manipulated to evolve atoms in a superposition of 
lattice site without giving time for the atom-atom inter- 
action, and the optical potential depth of the system can 
be configured just above the level where there is no direct 
tunneling. Then the dynamics of atoms in lattice will be 
dominated by the last term of the Hamiltonian in Eq. 

ignoring the atom-atom interaction. Therefore scal- 
ing up the scheme of using the QW on systems with large 
number of atoms in each lattice site or to infinitely large 
numbers of lattice site is also straight forward; whereas, 
the atom-atom interaction play a prominent role during 
the usual method of varying the potential depth. 



V. QUANTUM WALK WITH A NOISY 
CHANNEL AS A TOOLBOX 

Along with the coin parameters 9, ^ and a small 
amount of engineered noise can be used to control the 
redistribution of atom in the optical lattice. To demon- 
strate the effect of the QW with a noisy channel on atoms 
in optical lattice we consider a bit-flip channel, a phase- 
flip channel and an amplitude-damping channel. The 
bit-flip channel flips the state of the particle from |0) 
(|1)) to |1) (|0)) (Pauli X operation) and a phase-flip 
channel flips the phase of |1) to — 11) (Pauli Z opera- 
tion). We use the notation p for the noise level where, 
< p < 1. Therefore, the bit- (phase-)flip channel flips 
the state (phase) with probability (1 — p) during each 
step of the QW. An amplitude-damping channel leaves 
state |0) unchanged but reduces the amplitude of state 
|1) with probability (1 — p) resulting in an asymmetric 
distribution. The maximum decoherence effect using bit- 



and phase-flip channel is at p = 0.5 due to symmetries 
induced by these two channels during the QW; whereas, 
the amplitude-damping channel does not obey any sym- 
metry and hence the maximum effect is for p ~ 1 [l8|. 
Our numerical implementation of these channels evolves 
the density matrix employing the Kraus operator repre- 
sentation. 

Figure [7] is a comparison of the distribution obtained 
without noise channels to the distributions with phase- 
flip {p — 0.02 and p = 0.1) and amplitude-damping chan- 
nel {p = 0.2) on atoms initially in the MI state. Similarly, 
the redistribution of atoms in SF state when subjected 
to QW without and with noise channels is shown in Fig. 
in With the addition of phase flip noise of p = 0.02, a 
uniform distribution is obtained. Then the optical po- 
tential depth can be increased to cancel the correlation 
and obtain MI state. 

The effect of bit-flip on the distribution is close to the 
one obtained using phase-flip channel and hence the re- 
distribution in the Figs. [7] and [5] are shown only for 
phase-flip and amplitude-damping channels. Increasing 
noise level affects the variance of the QW 18] , which pro- 
portionally affect the atom-position correlation and the 
atom-atom overlap region. Therefore, it is important to 
restrict the noise level to a very low value {p ^ 0.1) to 
use noisy channel as a tool. 

Some of the distribution presented in Sees. IIVI and 
IVl using QWs are visibly similar to some of the distri- 
bution presented using different technique; in [s^l using 
time evolution density matrix renormalization group (t- 
DMRG) and in [s^, using quantum Monte Carlo sim- 
ulation. These technique present the change in density 
profile distribution with time, whereas in this paper we 
have discussed the change with number of steps of the 
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FIG. 8: Atoms in SF state after implementing a QW (A'' — 
M = 40,9 = 45°). With noiseless QW, the atoms spread, 
getting close to a uniform distribution. For phase damping 
Pp = 0.02 the distribution is uniform between ±20 and gets 
closer to Gaussian at pp — 0.1. Amplitude damping of state 
|1) introduces asymmetry to the distribution, pa = 0.2. 



nal state and the two states are transported in opposite 
directions. Delocalization of atoms over the seven lattice 
sites is demonstrated. This scheme can be adopted to 
implement the QW and to redistribute the density pro- 
file by introducing a r/-pulse ^] after each separation. 
An additional r/-pulse can be engineered to act as noise 
channels. 

In an implementation scheme proposed in Ref. [40| . 
a stimulated Raman kick, two selected levels of the atom 
are coupled to the two modes of counterpropagating laser 
beams to impart a translation of atoms in the position 
space. To realize transition from the SF to the MI phase 
using this scheme, atoms are first redistributed by imple- 
menting the QW in a long-Rayleigh-range optical trap 
and the optical lattice is switched on at the end to 
confine atoms and cancel the correlation. For transition 
from MI to SF, the atoms initially in the optical lattice 
are transferred to long-Rayleigh-range dipole trap and a 
QW is implemented. 



VII. CONCLUSION 



QW. 

VI. IMPLEMENTATION 

Optical lattices ranging from a simple periodic, square 
and cubic, to a more exotic ones, such as hexagonal and 
Kagome lattices using superlattice technique [37| have 
been created to trap and manipulate cold atoms. Ma- 
nipulation of cold atoms in time-varying optical lattice 
has also been reported [s^ . This has provided flexibility 
in designing and studying quantum phases and quantum 
phase transitions in coherent and strongly correlated ul- 
tracold atoms. 

In Ref. [3^ , controlled coherent transportation of ru- 
bidium atoms in the spin-dependent optical lattice has 
been experimentally demonstrated. Rubidium atoms in 
the MI state is evolved into superposition of two inter- 



In this paper, we have shown the use of the QW to 
study the dynamics of atoms in an optical lattice and 
expedite the process of quantum phase transition. We 
have also used the QW with experimentally realizable 
noisy channels to show the additional control one can 
have over the evolution and atomic density redistribu- 
tion. Theoretically, the evolution of the density profile 
with a QW can be used in place of quantum Monte Carlo 
simulation to study the correlation and redistribution of 
atoms in optical lattices. We expect the QW to play a 
wider role in simulating and expediting the dynamics in 
various physical systems. 
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